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ABSTRACT 



We present strong evidence of time stretching in the peak-to-peak time 
scales in the hght curves of BATSE Gamma Ray Bursts (GRBs). Extensive 
tests are performed on artificially dilated bursts to verify that the procedure 
for extracting the peak-to-peak time scales correctly recovers the stretching 
of bursts. The resulting robust algorithm is then apphed to the 4B GRB 
database. We derive a stretching factor of 1.92 ± 0.13 between the brightest 
burst group (P > 7 .7 photon ■ cm^"^ ■ s^^) and the dimmest burst group 
(P = 1.0 ~ lAphoton ■ cm~'^ ■ s~^) with several independent peak-to-peak 
time scale definitions and they agree within uncertainties. Such an agreement 
strongly supports the interpretation of the observed time stretching as 
time dilation caused by the cosmological expansion, rather than physical 
selection effects. We fit the result to cosmological models with A = 0, 
Qq from 0.2 to 1.0, and contraincd the standard candle luminosity to be 
Lq — 7.0 ± 2.0 ± 2.7 X lO^^photons ■ s~^. Our luminosity value is fully consistent 
with the value from the combined PVO and BATSE LogN-LogP curve with 
the BATSE bright bursts at low rcdshifts of Zbright = 0.11 ± 0.02 ± 0.025. This 
luminosity fit is definitely inconsistent with the the larger distance scale implied 
from associating burst density with star formation rates. 



Subject headings: gamma rays:bursts 
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1. Introduction 

Observations taken by the BATSE instrument aboard the Compton Gamma Ray 
Observertory have identified ~ 1800 GRBs and shown that their angular distribution is 
highly isotropic, and their distribution in space is inhomogeneous (Meegan et al. 1992, 
Briggs et al. 1996). Such a distribution naturally arises if GRBs are at cosmological 
distances. The logN-LogP distribution of GRBs is well studied and shown (e.g. Fenimore et 
al. 1993) to be consistent with cosmological models. The recent measurement of the redshift 
of GRB970508 (Metzger et al. 1997) between 0.835 ~ 2.1, and the possible identification of 
the host galaxy of GRB971214 (Kulkarni ct al. 1998) aX z 3.42 provide further evidence 
supporting the cosmological scenario. 

Piran (1992) and Paczyhski (1992) suggested that the light curves of the GRBs should 
be stretched due to the cosmological time dilation. This effect applies to all time scale in 
the GRB hghtcurves. In the past, different groups (Norris et al. 1994; Mitrofanov et al. 
1996; Bonnel et al. 1996; Rutledgc ct al. 1996 ) have investigated the correlation of the 
duration of bursts {T^q , Tgo, width of the main peak) and the burst brightness (peak fiux). 
The results have been contradictory, with different groups getting different results using 
identical method and similar data sets. 

Furthermore, questions have been raised whether the time stretching found is due to 
intrinsic correlation (Brainerd 1994, 1997) between pulse width and burst brightness for 
bursts drawn from a volume limited sample. It is argued that the observed correlation 
could arise either from the beaming of relativistic jets (Brainerd 1994), or the conservation 
of total energy in monoenergetic sources (Wijers & Paczyhski 1994). Thus the controversial 
claims of temporal stretching need neither imply cosmological distances nor be of any utility 
in understanding burst demographics. 
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In this Letter, we investigate the correlations of the time intervals between peaks with 
brightness indicators peak flux P. The peak-to-peak time scale is independent of earlier 
pulse width measurements. Attempts (Pozanenko et al. 1997, Norris et al. 1996) have 
been made to search for peaks separated by a valley with intensity difference of at least 
4(7. Such a definition is biased towards identifying more peaks in the bright bursts than 
in the dimmer bursts since the latter have less photon counts above background. We used 
artificially stretched bursts to identify a peak-to-peak time scale which correctly recovers 
the input stretching relations. 

2. Peak Finder Algorithm 

We have used the updated BATSE 64ms ascii database which contains 1252 bursts 
ranging from trigger 105 to trigger 5624. This database provides 64ms time bins from the 
concatenation of three standard BATSE datatypes, DISCLA, PREB and DISCSC. All 
three data types are derived from the on-board data stream of BATSE's eight Large Area 
Detectors (LADs), and all three data types have four energy channels, with approximate 
channel boundaries: 25 — bbkeV, 55 — llOkeV, 110 — 320keV, and > 320keV (see Fishman 
et al. 1989). The peak flux values are derived on the 256ms time scale as in the 4th BATSE 
catalog. 

The data were binned to 256ms to achieve better S/N ratio for the dim bursts. Bursts 
with Tgo shorter than 2s are excluded (Norris et al. 1994, Mitrofanov et al. 1996) from our 
analysis. The noise biases are rendered uniform (Norris et al. 1994) by diminishing the 
background subtracted signal to a canonical peak intensity of 2A00counts ■ s"^, and adding 
a canonical fiat background of SOOOcounts ■ s^^. The elimination of bursts with peak flux 
below the l.Ophoton ■ cm"^ • threshold is introduced to avoid trigger threshold effects 
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and ensure the statistical significance of the peaks identified. 

The detailed procedure to find all the peaks in a burst is as follows: (1) We fit the 
background using a quadratic function to the pre- and postburst regions. The background 
is subtracted from the data. (2) The maximum peak counting rate Cmax at 256ms is 
identified. (3) Identify all local maximums Cp that are separated by local minimums Cmin 
which satisfy Cp — Cmin > Ny ■ Cmax- (4) Each such local maximum also has to be greater 
than a threshold level, Cp > T ■ Cmax- is a parameter to ensure the peaks are distinct 
enough statistically. T" is a threshold level for accepting a local maximum as a significant 
peak. 

We used four definitions of peak-to-peak time scales. The first method uses all 
time intervals between successive peaks, Tj. Thus a burst with N peaks will provide 
N-1 peak-to-peak time scales, all of which have equal weight. Our second method is 
to logarithmically average all the time intervals between successive peaks Tp_p for each 
individual burst, resulting in one time scale per burst. Our third method is to consider only 
the time interval Ths between the highest and the second highest peak. The fourth method 
is to consider the time interval tfl between the first and the last peaks identified. 

3. Simulation Tests 

The intervals between peaks for long bursts (Tgo > 2s) varies from less than 1 second 
to tens of seconds, thus any dilation effect can be detected only in a statistical sense. 
Furthermore, the identification of peaks are complicated by small S/N ratios in weak bursts 
which might introduce noise peaks. The highly significant peaks do not suffer from these 
problems, while having more peaks identified will improve the statistics of our analysis. 
Therefore, it is not known a priori which parameters and T are the best choice to 
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measure the time intervals between peaks of GRBs, nor do we have a priori knowledge of 
which peak-to-peak time scale definition yields reliable results. A faithful procedure should 
be able to distinguish the time stretching effect from any systematic effects. 

To find a faithful procedure, we simulated the time stretching effects of the bursts and 
tested whether our procedure of using the peak finder algorithm in extracting the time scale 
correctly recovers the time stretching despite the systematic effects. Our procedure is to 
search the parameter space for and T to find regions where the peak finder algorithm 
recovers the input dilation with high confidence. The details of the simulation procedure go 
as follows: 

We assume a standard cosmology with A = 0, = "^Qo- Hence there is the following 
relationship between peak flux and redshift. 

^ 5.2xlO-5«-go^-(l + ^)-+^-L 

^' No + (go-i)(V2go^+i-i)F 

Here go is the deceleration parameter, a is the power law spectral index of GRBs, L is 
the burst luminsity in photons ■ s~^. In deriving this relationship, we assumed a Hubble 
constant of 75km ■ • Mpc~^. The power law spectral index varies from 1.0 to 3.0 for a 
majority of the bursts (Schaefer et al. 1994). We assumed go — 0.5, a = 2.0, and a typical 
(Horack et al. 1996, Hakkila et al. 1996) GRB power law luminosity function 4>{L) = AqL~^ 
{Lmin < L < L^ax, K = L^ax/Lmin = 100.0) with (3 = 2.0, Lmin = 2.3 X IQ^"^ photons ■ s~^ 
in the simulation process of stretched bursts. 

With this cosmology, we created a database of simulated BATSE light curves. Each 
burst was simulated to have peak flux value P randomly sampled from the LogN — LogP 
relation of the 4B catalog, and its luminosity L randomly sampled from the luminosity 
function 0(L) assumed. Both of them are sampled using the rejection method (Press et al. 
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1992). Its redshift z is then determined by Eq. 1. We then randomly select a burst from 
a library of 100 bright BATSE bursts, and determined its redshift z^ynght from the above 
relation. Burst photon counts from the 4 different channels from the Large Area Detectors 
are redistributed by a smooth interpolation function to redshift the photon energy. The 
energy shift factor is e = (1 + Zbright)/ (1 + -2)- 

For the burst at redshift z to be simulated, the relative time dilation factor is 
S = = (1 + z)/(l + Zfyright)- The light curve is then stretched by a factor of S. The 
photon count of ith bin Nj, in the stretched light curve is dependent upon the photon count 
of two bins Op and Op+i in the original light curve, where p is the integer portion oi i/S. 
The fraction of Op+i that will contribute is f = i — p ■ S {f = 1.0 ii i — p ■ S > 1.0). Taking 
account of the new background, the new count Ni should be: 

Ni = Op+i ■f + Op-{l-f)+ p{d), d = Op+i -f + Op-il- ff (2) 

where p(d) is a Poisson variable whose standard deviation is d^^^. The resulting light 
curve is then dimmed by a factor of P{z)/ P{zbright) as prescribed by Norris et al. 1994. 

This process is repeated 1250 times to generate a simulated database of gamma 
ray burst light curves, with the same LogN — LogP distribution as the real 4B GRB 
database. Each database simulation is repeated 10 times with different seeds for random 
number generation. By looking at the scatter between these 10 simulations we estimate the 
uncertainty associated with the simulation procedure. 

The peak finder algorithm is applied to each simulated database. We divide the 
bursts into 6 brightness bins with equal numbers of bursts, and logarithmically averaged 
the individual peak-to-peak times within each bin. The average peak-to-peak time scales 
derived are normalized to that of the brightest bin. The resulting relation is then compared 
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to the input dilation relationship of Eq. 1 by a test. The ^ values were calculated for a 
wide range of and T . 

In Fig. 1, the dependency on and T is shown by a contour plot for the simulation 
with the stretching relation of Eq. 1. We also performed simulations for bursts with no 
stretching, the dependency on and T is shown in Fig. 2. The contour areas that 
reside within x^ = 4 in each figure provide and T values that define a peak-to-peak 
time scale which faithfully recovers the input relation. Obviously, such much 
larger when we have no stretching (Fig. 2). The overlapping areas {Ny ~ 0.37, T < 0.22) 
provide a definition that reproduce the dilation relation with high confidence for both 
the stretched and unstretched simulated database. We repeated the simulation with a 
variety of cosmological parameters and luminosity values and found that such a conclusion 
consistently holds. In Table 1, the time stretching relation recovered from the algorithm 
is compared to the input relation. These simulation procedures ensure that we have a 
peak-finding algorithm which faithfully returns the stretching factors regardless of the size 
of stretching. 

The primary parameter that ensures the peaks found are nonstatistical is A^^. It 
is defined to be the peak-to-valley amplitude (measured in units of the highest peak 
height above background) separating each candidate peak from its neighbors. When iV^ 
is small, the algorithm picks up many statistical fluctuations as peaks, which distort the 
stretching relation. When is large, very few peaks are identified, and the statistics are 
poor. Nevertheless, these few peaks will result in large uncertainties if they are used as a 
measurement of such that Tj vs. P relation can be fitted to any model with a reasonable 
X^. In figures 1 and 2, this is evidenced by a region of x^ < 6 across the top for large 

values. While faithful, these regions carry little information as to distinguish between 
models. Therefore the algorithm only works for intermediate values. 
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The parameter T in the algorithm is defined to be a threshold level that all peaks 
identified have to exceed. When > T, the T constraint is primarily important to reject 
some large statistical fluctuations above and below the background. Although rare, such 
fluctuations usually yield a peak in the background data far away from the burst activity, 
creating large Tj values that are false and distort the apparent Tj vs. P relation. 

4. Results and Conclusions 

We applied the algorithm to the real BATSE database, and the result is shown in Fig. 3. 
A stretching factor S is calculated by fitting the resulting rvs.P relationship to Equation 1. 
The same procedure for simulation test is repeated for the other definitions of peak-to-peak 
time scales and the resulting algorithm is applied to the real BATSE bursts. The results are 
summarized in Table 2. The time stretching factor S of 1.92 ± 0.13 between the brightest 
(P > 7 .7 photon ■ cm~^ ■ s~^) and dimmest group (P = 1.0 ~ lAphoton ■ cm~^ ■ s~^) is 
found. Within the uncertainties, the different definitions of the peak-to-peak time scales 
give consistent results. 

Recently, two spectral class of bursts have been identified (Pendleton et al. 1997): those 
bursts with high energy emission (HE), and those with no high energy (NHE) emission. The 
latter group is predominantly fainter than HE bursts and is found to have a homogeneity to 
much lower fiux values. We applied the peak finder algorithm (see Table 3) to this subgroup 
with 232 bursts in the 3B catalog. Since the NHE bursts are inherently much fainter, there 
are fewer bursts with peak fiux P above the l.Ophoton ■ cm~'^ ■ s^^ threshold (66 out of 232) 
and hence the uncertainties are much larger. The time stretching factors of the NHE bursts 
appears to be consistent with those of the HE bursts (reduced xl = 1-01) and those of a 
no stretching model (reduced xl — 1-11). Therefore, it is safe to conclude that the limited 
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statistics available is not enough to determine whether any stretching exist in NHE bursts. 

It should be noted that the time dilation relation is independent of the GRB density 
evolution. We fit the result to the Eq. 1 and constrain the luminosity value Lq in a standard 
candle scenario. This provides an estimate of the redshift Zbright of brightest BATSE bursts 
independent of previous estimates from the LogN — LogP relation. We fit the data with 
varying values of deceleration parameters from 0.1 to 0.5, and power index of GRB spectra 
a from 1.0 to 2.0. We find the best fitting luminosity Lq = 7.0 ±2.0 ±2.7 x lO^^photons ■ s^^ 
(reduced xl — 1-23), and the brightest bursts (peak flux P — SOphotons • cm~'^s~^) at 
^bright — 0.11 ± 0.02 ± 0.025. The first error denotes the uncertainty associated with the 
varying values of qo and a, the second error is derived from the fits. This result favors 
the simple cosmological scenario of bright bursts at small redshifts rather than evolutionary 
scenario of bursts at much larger redshifts based on some theoretical arguments (Totani 
1997, Paczyhski 1997) that burst density traces star formation rates. 

It is well known that GRBs might have a broad luminosity function. We performed 
our simulation adopting a range dominated luminosity model (Horack 1996, Hakkila 1996) 
and demonstrated that for typical luminosity function (f){L) oc L~f^ {L^in < L < L^ax, 
K — Lmax/Lmin — 100.0) with powcr-law index /3 near 2, the broadening of the luminosity 
distribution will not smear the time dilation effect out. Our simulation shows that our 
technique can correctly recover the time stretching on a bright burst sample for this case. 
However, the time dilation data alone can not be used to constrain all the parameters 
Ljnin, Ljnax, P the powcr law luminosity function as it involves too many free parameters 
in the cosmological model. 

In summary, we have measured several peak-to-peak time scales of GRBs. Simulation 
tests are performed to verify that the time scales identified correctly reflect the physical 
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time scale intrinsic to the bursts. The ths a-nd '^fl time scales are similar to the durations, 
while the Tj and rp_p peak-to-peak time scales measurements are independent of previous 
duration and pulse width measurements. Our result unambiguously shows the existence of 
time stretching factor of s = 1.92 ±0.13 between the bright and dimmest GRB time profiles, 
and the agreement found between these multiple independent time scale measurements 
indicates that the temporal profiles of GRBs are universally stretched. Such an agreement 
implies that the time stretching relation comes from cosmological expansion rather than 
physical selection effects affecting particular property of bursts. 

We thank G. N. Pendleton for providing the hst of NHE bursts. This work has made 

use of the data obtained though the Compton Gamma Ray Observatory Science Support 
Center Online Service, provided by the Goddard Space Flight Center. 
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Pinhotons • cm, ^ • .s ^1 


^Nodilation 


^recovered 


1.0 ~ 1.4 


1.00 


1 04 ± 10 


1.4 ~ 1.8 


1.00 


1 05 ± 09 


1.8 ~ 2.5 


1.00 


1 13 ± 09 


2.5 ~ 3.6 


1.00 


0.98 ±0.08 


3.6 ~ 7.7 


1.00 


1.10 ±0.09 


> 7.7 


1.00 


1.00 ±0.09 


P{photons ■ cm~'^ ■ s~^) 


S Dilation 


^recovered 


1.0 ~ 1.4 


1.69 


1.82 ±0.12 


1.4 ~ 1.8 


1.56 


1.67 ±0.09 


1.8 ~ 2.5 


1.44 


1.58 ±0.11 


2.5 ~ 3.6 


1.30 


1.28 ±0.12 


3.6 ~ 7.7 


1.17 


1.12 ±0.11 


> 7.7 


1.00 


1.00 ±0.08 



Table 1: The input stretching factor and the measured factor using the definition of peak- 
to-pcak time scales < rp_p > with — 0.375 and T — 0.100 for both simulated stretched 
bursts and unstretched bursts. 



Peak-to-peak Timescale 
Definition 


Stretching Factor 

^observed 


n 


1.92 ±0.13 




1.96 ±0.17 


ths 


1.92 ±0.26 


tfl 


1.92 ±0.20 



Table 2: The observed time stretching factor between the brightest burst group (P > 
l.lphoton • cm~'^ • s~^) and the dimmest burst group (P = 1.0 ~ lAphoton • cm~'^ ■ s~^) 
using different definition of time scales. 
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P [photons ■ cm. ^ • s ^) 


snhe 


She 


sau 


1.0 - 1.4 


1.68 ±0.49 


1.83 ±0.23 


1.92 ±0.24 


1.4 ~ 1.8 


2.97 ± 1.28 


1.56 ±0.22 


1.65 ±0.20 


1.8 ~ 2.5 


0.90 ±0.28 


1.40 ±0.18 


1.27 ±0.16 


2.5 ~ 3.6 


1.54 ±0.60 


1.24 ±0.14 


1.42 ±0.16 


3.6 ~ 7.7 


1.40 ±0.70 


1.03 ±0.12 


1.08 ±0.13 


> 7.7 


1.00 ±0.29 


1.00 ±0.10 


1.00 ±0.09 



Table 3: The stretching factor relative to the brightest burst group for NHE bursts, HE 
bursts and all the bursts. Notice that the uncertainties associated with the NHE bursts are 
substantially larger due to the relative scarcity of such events in the 4B catalog. 
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Fig. 1. — The contour of the peak-to-peak algorithm for cosmological model with 
Lq = 2.3 X 10^'^photons ■ . There are two level of x^s, 4 and 6. The x and y axis 
are T and A^^, values respectively. The values are calculated by fitting rvs.P relation 
recovered by the algorithm to the input stretching relation of 6 burst brightness groups. 
There are 4 degrees of freedom for the fit, hence the region within — ^ provides a faithful 
algorithm which correctly recovers the input stretching relation. 
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>0.35 
0.3 
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Fig. 2. — The contour of the peak-to-peak algorithm for no dilation model. There are 
two level of x^s, 4 and 6. The x and y axis are T and values respectively. The values 
are calculated by fitting rvs.P relation recovered by the algorithm to the input stretching 
relation of 6 burst brightness groups. There are 4 degrees of freedom for the fit, hence the 
region within = 4 provides a faithful algorithm. Notice that the faithful areas are much 
larger in this case. 




Fig. 3. — The average peak-to-peak intervals of BATSE GRB database. Each dot in the 
graph represents the peak-to-peak interval identified Tj in a mutipeak burst. P is the peak 
flux in photons • cm~^ ■ s'^. Ti is the peak-to-peak interval in units of seconds. The solid curve 
is the best fit to a cosmological model with go — 0.2, A = 0, a = 1.5, and burst standard 
candle luminosity Lq = 7.0 ± 2.7 x lO^^photons ■ s"^. 



